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ABSTRACT 


An analytical study of catalyzed hydrazine decomposition reaction chambers 
was performed in order to establish procedures capable of predicting the steady- 
state and transient behavior of the system. Tne study included the development 
of computer programs to calculate the temperature and reactant concentration dis- 
tribution as functions of time and axial position in typical reaction chamber con- 
figurations. The computer programs evolved are based upon models of the reactor 
system which consider both thermal and catalytic decomposition of reactants, along 
with simultaneous heat and mass transfer between the free -gas phase and the gas 
within the pores of the catalyst pellets. Results obtained using these computer 
programs are. in excellent agreement with data available from small scale engine 
firings using Shell 405 catalyst thus permitting not only correlation of the test 
data but also indirect determination of hydrazine and ammonia decomposition kine- 
tics on 405 catalyst. 
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FOREWORD 


This work was performed by United Aircraft Research Laboratories for the 
National Aeronautics and Space Administration under Contract NAS 7-^58 initiated 
April 15., 1966 . 

Included among those who cooperated in performance of the work under Contract 
NAS 7-^58 were Dr. A. S. Kesten, Program Manager, Dr. W. G. Burwell, Chief, Kine- 
tics and Thermal Sciences Section, Mr. D. B. Smith, and Mrs. E. Smith of UARL. 

This work was conducted under program management of the NASA Chief, Liquid 
Propulsion Experimental Engineering Systems, NASA Headquarters, Washington, D. C., 
and the Technical Manager was Mr. T. W. Price, Jet Propulsion Laboratory, Pasadena, 
California 
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Analytical Study of Catalytic Reactors 
for Hydrazine Decomposition 
First Annual Progress Report 
April 15, 19 66 - April 14-, 1967 
Contract No . MAS 


SUMMARY 


The Research Laboratories of United Aircraft Corporation under Contract NAS 
7-4-58 with the National Aeronautics and Space Administration are performing an 
analytical study of catalytic reactors for hydrazine decomposition. This first 
annual technical report summarizes work performed under this continuing contract 
from April 15., 1966 to April 14- , 1967 * Work during this period has included the 
development of computer programs which are used to calculate the temperature and 
concentrations of hydrazine and ammonia as functions of time and axial position in 
typical reaction chamber configurations. The computer programs are based upon 
steady-state and transient models of the reactor system which consider both thermal 
and catalytic decomposition of reactants, along with simultaneous heat and mass 
transfer between the free-gas phase and the gas within the pores of the catalyst 
pellets . The steady-state model contains a detailed mathematical description of 
the processes of heat and mass transfer in the pores of the catalyst pellets, while 
the transient model uses over-all coefficients to evaluate the effects of heat and 
mass diffusion on reaction rates in the porous catalyst particles. These coeffi- 
cients are defined in terms of the physical properties and operating parameters of 
the reactor system. 

The computer programs representing the steady-state and transient models of 
the reactor system have been used to calculate temperature and reactant concen- 
tration distributions for various initial bed temperatures, feed temperatures, 
chamber pressures, mass flow rates, injector locations, and catalyst size distri- 
butions. The results have been compared, wherever possible, to available measure- 
ments of temperature profiles and exhaust concentrations made during firings of 
typical hydrazine engine reaction chambers employing Shell 4-05 catalyst . Excellent 
agreement between theoretical and experimental results is indicated. 

This agreement has permitted verification of kinetics information for use in 
both the steady-state and transient programs. In particular, estimates have been 
permitted of the rate laws and associated rate constants governing the heterogen- 
eous, vapor-phase decomposition of hydrazine and ammonia on the Shell 4-05 catalyst. 
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UHRODUCTION 


Effective design of distributed-feed catalyzed monopropellant hydrazine rocket 
engines and gas generators requires accurate procedures for predicting the effects 
of the design parameters of the reactor system on the steady-state and transient 
performance of the system. Prior to the present undertaking, this general capa- 
bility did not exist, although simplified methods had been developed for predicting 
the steady-state performance of idealized catalytic reactors (Refs. 1 and 2). These 
simplified methods do not adequately describe the combined processes of heat trans- 
fer, diffusion, and chemical reaction and are restricted to consideration of a sin- 
gle decomposition reaction at any axial location within the reaction chamber. To 
overcome these deficiencies, a more comprehensive theoretical analysis was deemed 
necessary. Preliminary investigations at United Aircraft Research Laboratories 
demonstrated the feasibility of such an analysis and, further, indicated that both 
transient and steady-state performance characteristics could be predicted. 

Based upon the above investigations, a comprehensive analytical program was 
formulated with the objectives of (a) developing computer programs for predicting 
the temperature and concentration distributions in monopropellant hydrazine cata- 
lytic reactors, and (b) performing calculations using these computer programs to 
demonstrate the effects of various system parameters on the performance of the 
reactor and, ultimately, to establish guidelines for reactor design. During the 
first year of the contract period attention has been focused on describing the 
steady-state and transient behavior of a reactor system in which complete radial 
mixing in the free-gas (or liquid) phase was assumed. Currently, the steady-state 
program is being extended to include radial as well as axial variations in temper- 
ature and concentrations. This portion of the program, scheduled for completion 
in April of 1968, will permit an analysis of various injection schemes which exhibit 
radial nonuniformities . Effort during the first annual reporting period is des- 
cribed in detail in succeeding sections of this report. 
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DISCUSSION 

The analysis of a hydrazine engine reaction system carried out to date per- 
tains to a reaction chamber of arbitrary cross section packed with catalyst par- 
ticles into which liquid hydrazine is injected at arbitrarily selected axial 
locations. At these locations, hydrazine injection is taken as uniform across the 
cross section of the chamber. Catalyst particles are represented as "equivalent" 
spheres with a diameter taken as a function of the particle size and shape. Both 
thermal and catalytic vapor phase decomposition of hydrazine and ammonia are con- 
sidered in developing equations describing the concentration distributions of these 
reactants. Diffusion of reactants from the free-gas phase to the outside surface 
of the catalyst pellets is taken into account. Since the catalyst material is 
impregnated on the interior and exterior surfaces of porous particles, the diffusion 
of reactants into the porous structure must also be considered. In addition, the 
conduction of heat within the porous particles must be taken into account since the 
decomposition reactions are accompanied by the evolution or absorption of heat. 

A treatment of the transport processes described above constitutes a general 
model of the reaction chamber. The analysis of the steady-state performance of a 
reactor based on this model can be carried out in a straight-forward manner. How- 
ever, the analysis of the transient behavior of the system using this general 
model is quite complicated since both the reactant concentrations and the rate 
constants for decomposition, which are exponentially dependent on temperature, will 
be functions not only of time, but also of position within the catalyst particles. 

To circumvent these complications, the transient behavior of the system can 
be analyzed by employing a second model which takes heat and mass diffusion within 
the catalyst particles into account by considering their effects on the rates of 
the catalytic reactions through the "utilisation factor." This factor, which is 
defined as the ratio of actual reaction rate to that which would be produced by 
the catalyst particle if negligible temperature and concentration gradients 
existed within the particle, can be calculated as a function of temperature and 
reactant concentration at the outside surface of a catalyst particle at steady- 
state. If it is assumed that, even under transient conditions, utilization factors 
are functions only of surface temperatures and concentrations, the analysis of the 
transient behavior of the reactor system becomes tractable. Using this method then, 
the temperature and reactant concentrations can be described as functions of time 
and position in the reaction chamber. 

Included in succeeding sections of this report are detailed descriptions of 
(a) the development of the computer programs representing the steady-state and 
transient models of the reactor system, (b) the establishment of kinetics infor- 
mation pertinent to Shell U 05 catalyst particles (Refs. 3 , 5, and 6) for use in 

both of these programs, and (c) the use of these programs to calculate temperature 
and reactant concentration distributions for various reactor operating conditions. 
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Steady-State Model 

In developing the present steady-state model, the temperature and reactant 
concentrations in the interstitial phase (i.e., the free-fluid phase as distin- 
guished from the gas phase within the porous particles) are assumed to vary only 
with axial distance along the bed. In the entrance region of the reaction chamber, 
where the temperature is low enough to permit the existence of liquid hydrazine, 
vaporization of liquid is assumed to occur as a result of decomposition of vapor 
hydrazine within the pores of the catalyst particles. That is, catalytic reaction 
is assumed to be fast enough to keep liquid hydrazine from wetting the pores of 
the particles; the hydrazine concentration at the surface of the catalyst particles 
at any axial location in the entrance region is then computed from the vapor pres- 
sure of liquid hydrazine in the interstitial phase at the same axial location. 
Neglecting axial diffusion of heat or mass, the change in enthalpy of the inter- 
stitial phase in the region where liquid hydrazine is present (i.e., where hj < h, ) 
is related to the concentration gradient at the surface of the porous catalyst 
particles by 


G 


d hj 
dz 


+ 


H N 2 H «d a 

n u p M p 


(dc^ HA \ 


+ 


F (h. - hj 

I F 


0 for h| < h| V (l) 


The variation of mass flow rate, G , with axial distance is easily computed from 
the rate of feed of liquid hydrazine from the distributed injectors into the system. 
In the region where liquid hydrazine exists at temperatures below the vaporization 
temperature, the temperature may be obtained from 


T i = T f + 


h i - h F 


'f 


for hj < hj 


( 2 ) 


In the two-phase region, where Tj = T vap , the weight -fraction of vapor may be com- 
muted from 


WEIGHT - FRACTION VAPOR 


hj - hj 

Tv 

h, - hj 


for hj L < hj < hj V 


(3) 


At the axial position at which the enthalpy of the interstitial phase is just equal 
to the enthalpy of vapor hydrazine at the boiling point ( hj - hj V ), the fraction 
of hydrazine injected upstream of that point which has been decomposed is easily 
calculated from an overall heat balance. The associated amounts of ammonia, nitro- 
gen, and hydrazine formed from decomposition of hydrazine can then be calculated 
taking the decomposition reaction as 


N 2 H 4 


2NH 3 + N 2 4- H 2 


4 


F910461-12 


It should be noted t ha t this is the overall reaction scheme determined experimen- 
tally for both homogeneous decomposition of hydrazine (Refs. J, 8, 9) and low 
pressure heterogeneous decomposition of hydrazine on platinum surface 


(Rof mV 

V'"' 4 - • —T / * 


In the remainder of the reaction chamber, where hj>hj V , heat is being sup- 
plied to the system by homogeneous as well as heterogeneous decomposition of hydra- 
zine . In addition, at sufficiently high temperature, heat is removed from the 
system by the endothermic decomposition of ammonia. For h j > h j v then, the change 
in enthalpy with axial distance is related to the concentration gradients at the 
surface of the porous catalyst particles by 


dz 


n 2 h 4 

+ H r. y 

hom 


+ A r 



N 2 H 4 


+ 



F ( h | - h p ) 


0 


( 4 ) 


The changes in reactant weight fractions in the interstitial phase with axial 
distance are related to 'the concentration gradients at the surface of the porous 
catalyst particles by 


dz 



r hom^ 





FCj N a H « 

G P i 


( 5 ) 


^nh 3 

dz 



(6) 


*Tt is more commonly assumed, without benefit of experimental evidence, that the 
decomposition reaction is 3 Np^4 — ^ ^ NH3. + Rjj* followed by dissociation of one of 
the four ammonia molecules to nitrogen and hydrogen. This two-step process leads 
to the same overall reaction cited above but assumes that a minimum of 25 percent 
of the ammonia produced by hydrazine decomposition also decomposes. The fractional 
ammonia dissociation, f , calculated assuming the validity of the two-step process 
is related to the fractional ammonia dissociation calculated in the present report 

by 


( f ) two-step = 3(f) present + 1 

process report 

4 


5 


F9io46i-i2 



The changes in reactant concentrations with axial distance are then given by 


otj J d Wj J C| J dp\ 

dz dz + pi dz 

where 

\ dpi \ du i dii 

P\ dz M dz Tj dz 

and 


(9) 


(ic) 


j_ dM 
m dz 




dw i_J 
dz 


(ll. 


The temperature of the interstitial phase in this region is related to the 
enthalpy by 


h i 



/ T ' c f dJ, < 

T 

vop 


( 12 ) 


where Cf is the specific heat of the gas mixture and is a function of temperature 
as well as the concentration of the constituents of the mixture . 


The reactant concentration profile in the porous particles at any axial loca- 
tion must satisfy the diffusion equation for mass transport as well as an analogous 
equation for heat conduction. Neglecting the effect of the translational motion 
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of the gas stream on diffusion within the porous particles and assuming constant 
diffusion coefficients, D p , and thermal conductivities, K p , these equations may be 
written as 


D_ V 2 C. - r 


het 


0 


(13) 


Kn v 2 T n - Hr 


het 


0 


(14) 


Using Eqs. (l3) and (l4), Prater (Ref. 11) has pointed out that temperature and 
concentration are related quite simply by 


( t pI 


HD p 

Kn 


(Cp) s - Cp 


(15) 


The use of this relationship enables the reaction rate, r het , to be written as a 
function of concentration alone instead of concentration and temperature . Equation 
(13) can then be solved for the concentration at any point in the porous particle 
in terms of the concentration at the surface of the particle, (C p ) s . The solution 
is derived in Appendix A as an implicit integral equation given by 


Cp( x) r (C p ) s 




r het ( c p ) wr 
D p 




eZ r het ( c p) 
D P 


dt 


(ie 


The rate of chemical reaction on the catalyst surfaces, r het (C p ) , is given, in 
general, by (Ref. 12) 


r het (c p ) - k 0 (c p )' s n (c p ) n exp {x/3 (I - c p /(c p ) s )/[ I + (I - c p / (c p ) s ) ]} (17) 


where the parameters, k 0 , n , / , and /3 are defined in the list of symbols. 

If the reactant concentration at the surface of the particle, (C p ) s , is equal 
to the concentration in the interstitial phase, Cj , Eqs. (l6) and (17) can be 
solved simultaneously with the equations describing the changes in enthalpy and 
reactant concentrations in the interstitial phase with axial distance to yield the 
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steady-state temperature and concentration profiles in the reaction chamber. In 
the liquid region of the chamber, for example, it is assumed that liquid hydrazine 
wets the outside surface of the catalyst particles so that (Cp) s NzH4 : C \ NzH4 

where c ; N ? H4 is the vapor concentration in equilibrium with liquid hydrazine at 
temperature Tj . In the vapor region of the reactor it can be shown that (C p ) s 3 ~ 
c j nh 3 at the temperatures existing in the chamber since the rate of transport of 
material from the bulk fluid to the outside surface of the catalyst particles is 
much greater than the rate of dissociation of ammonia within the particles. How- 
ever, in the case of hydrazine decomposition in the vapor region the rate of cata- 
lytic decomposition is so high that (Cp)^ 2 ^ « Cj NzH4 and the reaction rate is 
governed by the rate of mass transfer from the bulk vapor to the outside surface 
of the catalyst particles. This rate is given approximately by (k c Cj) 2 4 where 
the mass transfer coefficient, k c N2H4 , maybe estimated from (Ref. 13) 



0.61 G 

P\ 



G 

A p 



(18) 


[ I n 2 h 4 

For hydrazine decomposition in the vapor region, then, the terms D p (o'Cp/o'x) 
in Eqs . (4), ( 5 ), (6), ( 7 ) and (8) can be replaced by (k c Cj) N 2 H 4. 

In the liquid-vapor region the situation is somewhat more complicated since 
it is difficult to predict whether liquid or a combination of liquid and vapor wets 
the outside surface of the catalyst particles. Both of these options are presently 
in the computer program representing the steady-state model. In the case in which 
both the liquid and vapor are taken to wet the particle surface, it is assumed that, 
at a given axial location, the fraction of the surface covered by vapor is equal 
to the weight -fraction of vapor present. Decomposition rates, computed assuming 
pure liquid surface coverage and then pure vapor coverage, are weighted accordingly. 
Fortuitously, for the system considered here, the liquid-vapor region is so narrow 
that the choice of either of these options has negligible effect on the resulting 
temperature distributions. This point is discussed further in the section covering 
the results of steady-state calculations. 

Finite difference methods have been used to program for digital computation 
the ordinary differential equations describing the changes in enthalpy and reactant 
concentrations in the interstitial phase. No iteration is necessary to solve these 
equations numerically when the incremental axial distances are sufficiently small. 
The size of a succeeding increment is calculated at each axial position as a func- 
tion of the rates of change of temperature and fractional ammonia dissociation with 
axial distance. However, Eq. (l6), which must be solved simultaneously with the 
differential equations, is an implicit integral equation which requires an itera- 
tive procedure for solution. Hand calculations have indicated that convergence to 
a solution for C p (x) is difficult to achieve unless the initial estimate of the 
concentration distribution is fairly accurate. Methods have been developed for 
generating this estimate and an iterative procedure has been devised which effects 
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rapid convergence over a fairly wide range of conditions . This procedure is pre- 
sently used as a subroutine in the main program representing the steady-state model. 
Sample calculations have been made to illustrate the use of this subroutine for 
computing concentration profiles in the porous catalyst particles, the associated 
temperature profiles, and the concentration gradients at the particle surface. The 
subroutine was run for values of the reaction rate constant, k 0 > between 2.0 sec -1 
and 2.0 x 10° sec" 1 , and for (3 = 0.171 } Y = 18.3, Dp = 0.84 x 10“ 5 ft 2 / sec, a = 
0.005 ft, n = 1, and (c p ) s = 0.357 Ib/ft^. The values chosen for fi , y , Dp, 
n , and (Cp) s are reasonable for the decomposition of hydrazine within Shell 405 
catalyst pellets . The resulting concentration profiles are plotted for each reac- 
tion rate constant in Fig. 1 and the associated temperature profiles are shown in 
Fig. 2. The concentration gradient at the particle surface is plotted as a func- 
tion of the reaction rate constant in Fig. 3» 

Transient Model 


In developing the transient model, the temperature and the concentrations of 
reactants in the interstitial phase are assumed to vary only with time and axial 
distance along the bed. In this system film coefficients are used to describe heat 
and mass transfer between the interstitial phase and the outside surface of the 

catalyst pellets. The reactant concentrations, c D , and the temperature, T p , are 
taken as uniform within the interior’ of the porous particles. Heat and mass dif- 
fusion within the particles are taken into account by considering their effects on 
the utilization factor previously defined. That is, reaction rates computed on the 
basis of uniform c p and T p are multiplied by a utilization factor determined by 
analogy with the steady-state system. In addition, it is assumed that liquid velo- 
cities are sufficiently low relative to other rate processes so that, for all prac- 
tical purposes, steady-state in the liquid and liquid-vapor regions is achieved as 
soon as the liquid reaches a given axial location in the reactor. 

The transient model is concerned then with the vapor region only, where velo- 
cities are about three orders of magnitude greater than in the liquid region. Here 
it is reasonable to assume that gas velocities are so great that the time lag from 
the entrance to the vapor region to any position z for the fluid is negligible com- 
pared with other transient effects. With this assumption the rates of change of 
temperature and reactant weight fractions with axial distance in the interstitial 
phase are given by 

dlj_ _ _ H n z h « g _ _A £_ 
dz ' GCp rhom GC f 

(hj-hp) + Cp (Tj ~ivap) 


F 

GC F 




T p) 


1 1 
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j_ j 

dz 6 I 


hom 


8 — An k 


p 


N 2 H 4 


( c 


NoH 


2 n 4 _ 


N?H 


2"4 


- N_H. 

hC; - 


G P i 


( 20 ) 


aw, NH ^ 

dz 


G 



M 


nh 3 


M 


n 2 h 4 


. . NH 3 

Apkc 


( c 


NH, 



F C| NH 3 
G Pi 


( 21 ) 


n 2 

aw, z 

dz 


l_ ( s M Nz 
G 1 rhom 2M^ h ^ 


An kc 



C 



F C j N 2 
G pi 


( 22 ) 


aw , H2 I I M h z 

dz ' G | r hom ° 2 M N 2 H4 



F C j H 2 

G Pi 


(23) 


The changes in reactant concentrations with axial distance can then be obtained 
using the above equations together with Eqs. (9)> (lO), and (ll). At a given 
axial location the rates of change of temperature and reactant concentrations in 
the catalyst particles with time are given by 


tfT p 

d\ 


i 


Ps 





a Ps c s 


( T i - t p) 


(24) 


d c 


N 2 h 4 

P 


d t 



n 2 h 4 

het 


3k c N 2 H4 


a 0 a 



Cp N * H *) 


(25) 


d t 


N->H 


2 n 4 


het 


M™* 

M* 2 " 4 


r nh 3 

het 


+ 




(26) 
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dc p Nz 
d\ 


a P 


N 2 H 4 

r het 2M N 2 H 4 


NH 3 M N 2 

r het 2 m nih 3 


3k c Nz 
a p a 




(• 27 ) 


dc " 2 

d\ 



n 2 h 4 

net 


M h 2 


2M 


n 2 h 4 


+ 


r nh 3 

r het 


3M h 2 n 



3k c Hg 

a p a 



(28) 


Recalling that the reaction of hydrazine on the catalyst surfaces is extremely fast, 
so that the reaction rate is controlled by the rate of transport of hydrazine to 
the catalyst surfaces, Eq. (25) can be used to define r het N 2 H 4 by no ting that 
(<7Cp N2H4 /<7t ) and C p N 2 H4 are both approximately equal to zero. The reaction rate is 
then given by 


3k 


n 2 h 4 


het 


(29) 


The reaction rate of ammonia on the catalyst surfaces, r het NH3 , is computed by 
multiplying the rate of reaction calculated on the basis of uniform T p and c p 
(see section on Kinetics Information) by the utilization factor determined by 
analogy with the steady-state system. Under steady-state conditions the concen- 
tration gradient at the particle surface, ( c/Cp ^ 3 / dx ) , can be used to define 
this utilization factor which, for a first-order reaction, is given by: 


UTILIZATION FACTOR 


3 D P (<7c P /c/x)] c 


( c p)gko] 


(30) 


The subroutine in the steady-state program which calculates concentration distri- 
butions in the catalyst particles has been used to compute the concentration gra- 
dients at the particle surface and the corresponding utilization factors for about 
100 different reactor operating conditions covering the range of interest of each 
of the parameters in this program. The result s were then plotted as a function of 
a dimensionless diffusion parameter, a Vko/D p . All of the results fell within the 
limits of the dashed curves shown in Fig. 4. The mean results are shown as the 
solid curve in this figure. 

Heat and mass transfer coefficients may be estimated from (Ref. 13 ). 


0.61 G \ / jx \-0 667 / G \-0.4l 
Pi / Ui Dj J / Up^J 


(31) 
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(32) 


The ordinary and partial differential equations representing the transient 
model of the reactor system (Eqs. (19) through (28) along with Eqs. (o'i through 
(ll)) must be solved simultaneously for the variation of temperatures and reactant 
concentrations with time and axial position in the reactor. The partial differ- 
ential equations representing the changes in interstitial temperature and reactant 
concentrations with axial distance can be treated as ordinary differential equa- 
tions by integrating them at constant time. This can be accomplished, together 
with integrating Eqs. (2 h) through (28) at fixed axial positions, by establishing 
a network containing fixed time intervals and intervals of axial distance . Ideally, 
then, reactant concentrations and temperatures in the interstitial phase and in the 
catalyst particles could be computed at the 

(U . = (Tp). . + 

',1 i - 1 . ) 


1 th time and j th axial distance from 


dir. 


d t / i-i, j 


At 


(33) 





(34) 





H 


Az 


(35) 



(36) 


This procedure is satisfactory as long as the intervals of time and axial distance 
are chosen to be sufficiently small so that the differential quantities are approx- 
imately constant in each interval. In the system considered here this restriction 
leads to rather long computing times. The restriction can be relaxed somewhat by 
rearranging each of the differential equations in the form 


-^ = a - £g (37) 

ds 
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where it is required now that only the quantities a and /3 remain approximately 
constant while integrating the equation from S k _ , to S k (corresponding to At or 
Az )• Equation (37) can be integrated to obtain 


9k 



s k-i) a 
+ £ 


g P (S k “ S|(.| ) 


(38) 


where g k is the value of g at s k , and g k . ( is the value of g at s k .| • An 
alternative form of Eq. (38) is 


’ k-i 


+ /^a\ 

r - /3 ( S k - S k _,) 

1 - e 

1 Ws/k-i 

£ 


(39) 


It can be seen that, for small values of ( s k - S k _ ( ), Eq. ( 39 ) reduces to the same 
form as Eqs. ( 33 ) through ( 36 ). It is convenient to use equations of the form of 
Eq. (38) in place of Eqs. (33) and ( 3 A) to compute particle concentrations and 
temperatures, and to use equations of the form of Eq. ( 39 ) in place of Eqs. (35) 
and (36) to compute interstitial concentrations and temperatures. 


Kinetics Information 


A number of investigators have reported rate expressions for the kinetics 
of homogeneous, vapor-phase decomposition of hydrazine (Refs. 7, 8, and 9)- The 
rate expression chosen for use in this work was that of Ref. 7 because the results 
were obtained in a reactor similar to the one of interest here. The equation 
reported is 


N 2 H 4 

r hom ' 


2.14 x lO 10 


c . n 2 H 4 q - 33,000/ Tj 


lb/ft 3 - sec 


(Ao) 


where Cj N 2 H4 is in lb/ft/ and Tj is in degR. 

Similar information relating to the catalytic decomposition of hydrazine is 
not available. However, a rate expression of the same form as Eq. (Ao) was chosen 
to represent the kinetics of the heterogeneous decomposition. Qualitative evidence 
from rocket firings at various feed temperatures indicates a rather low activation 
energy. A value of 2500 x 1.986 Btu/lb-mole was chosen rather arbitrarily and a 
preexponential factor in the Arrhenius type rate law was computed by using engine 
test data to make a rough estimate of the reaction rate at a typical feed temper- 
ature. These approximations have led to the following expression for reaction 
rate: 


r het = l0 '° (? ~ 250° / t p 


lb/ft 3 - sec 


(Al) 


where Cp^ 4 is in lb/f/ and T p is in deg R. 
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Melton 
experiments 
They report 


(Ref. i4) and Logan and Kemball (Ref. 15) have conducted low pressure 
on the vapor-phase decomposition of ammonia on platinum catalysts, 
a reaction rate expression of the form: 


nh 3 

nh 3 _ nh 3 c P - 50,000 /T P 

r het ' a H 2 n e 

[Cp 


lb/ft 3 - sec 


( 42 ) 


where T p is in deg R and n is given as 1.0 by Melton and as 1.6 by Logan and 
Kemball. Melton states however that "... The inhibitive action of hydrogen is 
probably a little greater than that corresponding to the first power of hydro- 
gen. . ." A similar dependence on hydrogen concentration at higher pressures was 
questionable particularly if hydrogen was physically absorbed on the platinum cata- 
lyst. Increasing hydrogen concentration beyond a certain point might then have 
much less effect on impeding the reaction rate. It was decided to perform computer 
calculations using the steady-state program assuming no dependence of r het 3 on 
the hydrogen concentration and then see if the value of a NH3 which best fits the 
experimental data decreases as the pressure increases. Such behavior would indi- 
cate that the dependence on hydrogen concentration did indeed extend to higher 
pressures. This behavior was in fact observed. It was assumed, of course, that the 
order of the ammonia dissociation reaction with respect to hydrogen is the same oh 
the Shell 405 catalyst as it is on platinum. Although this assumption remains un- 
tested, subsequent calculations using the steady-state program exhibit the best 
agreement with the available engine test data when r ^ e j NH3 : '- s §-^ ven ^ 


r het = °' 3x| ° 


r nh 3 
C P 

1.6 


(c P Hz ) 


e 


30,000/ Tp 


lb/ft' 


sec 


( 43 ) 


where the concentrations are expressed in lb/ft^. 

It should be emphasized that the value of 0.3 x 10 11 for a^ 3 was obtained 
by comparing experimental data with temperature and concentration profiles calcu- 
lated using the computer program representing the steady-state model. It is impli- 
citly assumed in this model that the temperature at the surface of the catalyst 
pellet is equal to the interstitial temperature at the same axial location. This 
assumption has little effect on the calculated rates of the diffusion-controlled 
hydrazine decomposition and it is fairly accurate in regions of the reaction chamber 
where only ammonia is decomposing. However, in regions where both hydrazine and 
ammonia are decomposing simultaneously, the catalyst particle surface temperatures 
may be higher than the interstitial temperatures . If this were taken into account 
in the steady-state model the value of a NHa exhibiting closest agreement with 
experimental data would be somewhat larger than 0.3 x 10 11 . In the simpler format 
of the transient model the differences between the particle surface and interstitial 
temperatures are accounted for and the rates of heat transfer corresponding to these 
temperature differences are expressed in terms of a heat transfer coefficient. 
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Theoretically then, the transient model, when extrapolated lo steady-state and 
compared with experimental data, should yield a more accurate value of q NH 3 . The 
value obtained using the extrapolated transient results is 1 x 10 -*-l. However, 
heat transfer coefficients used in this model are estimated using correlations devel- 
oped for nonreacting systems . In regions where significant chemical reaction takes 
place , such as the zones in which both hydrazine and ammonia are decomposing, actual 
heat transfer rates are considerably higher than those calculated in the transient 
model. This would lead to a value of a NHj less than 1 x 10 11 . The actual value 
of a 3 is probably between the one estimated using the steady-state program 
(0-3 x lCp-l) and the one estimated using the transient program (l x 10^-4). 

Results of Calcinations 


Calculations were made with the computer programs representing the steady- 
state and transient models of the reactor system to evaluate the effects on temper- 
ature and reactant concentration distributions of such parameters as initial bed 
temperature, feed temperature, chamber pressure, mass flow rate, axial injection 
profile, and catalyst size distribution. The sensitivity of the calculated results 
to variations in the rate constants associated with the catalytic decomposition 
reactions was also determined. The results have been compared, wherever possible, 
to available measurements of temperature profiles and exhaust concentrations made 
during firings of typical hydrazine engine reaction chambers. 

Steady-State Results 


A series of calculations was made using the steady-state program for cases 
for which experimental information is available (Ref. 2 ) in order to examine the 
effectiveness of the steady-state model and to establish the rate constants associa- 
ted with the catalytic decomposition reactions. The results of these calculations 
are shown in Figs. 5 through 19- For each of these cases the catalyst bed packing 
was taken to consist of 25-30 mesh catalyst particles for the first 0.2 in. and l/8 
in. x 1/8 in. cylindrical pellets for the remainder of the bed. This configuration 
is referred to in the figures as the "standard bed configuration." Except where 
the rate constants a NH3 and a N2hl4 were varied, the rates of the catalytic decom- 
position reactors were given by Eqs . (4l) and (4-3). Equation (43) was simplified 
somewhat by noting that, since hydrogen diffusion is quite rapid, Cp Hz ~ C j Hz . 

The reactor operating conditions of chamber pressure and mass flow rate were varied, 
as was the feed temperature, T F ; although, in each of these cases, T F was between 
approximately 510 and 530 deg R. 


Temperature distributions are plotted in Fig. 5 for a case in which the reactor 
operating conditions were taken as G = 3-12 lb/Tt^- sec and P = 479-5 psia. Temper- 
ature profiles are plotted for a Nhl3 = 0.3 x 10l0, 0.3 x lcA-1, and 0.3 x 1CA-2 
(ib/ft 3 ) °'^ (sec) -1 . An enlarged view of the temperature profile in the liquid- 
vapor region is shown to illustrate the effects of pure liquid as opposed to a 
combination of liquid and vapor wetting the outside surface of the catalyst parti- 
cles. Also shown in Fig. 5 are temperature measurements obtained by Rocket Research 
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Corporation (Ref. 2) during the course of engine firings under the same operating 
conditions. Comparison of the calculated results with the experimental data indi- 
cates better correspondence of results with a NH3 = 0.3 x 10-^ (lb/ft3)0-6 (sec) - ^. 
The calculated mole-fraction profiles associated with this rate constant, together 
with experimental values of mole -fract ions , are shown in Fig. 6. The corresponding 
fractional ammonia dissociation profile is shown in Fig. J. The results of similar 
calculations made for G = 2.43 lb/ft 2 -sec and P = 1042 psia are shown in Figs. 8 
through 10; those made for G = 6.29 lb/ft 2 -sec and P = 974 psia are shown in Figs. 

11 through 13; those made for G = 1.52 lb/ft 2 -sec and P = 217-9 psia are shown in 
Figs. 14 through l6, and those made for G = 1.51 lb/ft 2 -sec and P = 111.4 psia are 
shown in Figs. 17 through 19. 

In each of the above calculations the rate constant for the catalytic decom- 
position of hydrazine was taken as 10 10 sec -1 . Since catalytic decomposition of 
hydrazine is diffusion-controlled in the vapor region, varying this rate constant 
would only affect the width of the liquid and liquid-vapor regions . These regions 
are extremely narrow for the cases examined above because of the large surface area 
of the small catalyst particles in these regions and because of the high vapor 
pressure of liquid hydrazine at the feed temperatures considered. 

The effects of various reactor operating conditions on steady-state temper- 
ature profiles are illustrated more specifically in Figs. 20 through 23. A 
reference case was chosen in which the operating conditions were G = 3.0 lb/ft 2 - sec, 
P = 100 psia, and T F = 530 R, and the catalyst bed configuration was the 

"standard bed configuration." Chamber pressure, mass flow rate, catalyst bed con- 
figuration and axial injection profile were then varied in turn and the resulting 
temperature distributions were plotted. In Fig. 20 temperature distributions are 
plotted for chamber pressures of 100, 500, and 1000 psia with all other conditions 
taken as those of the reference case. Increasing pressure causes the peak temper** 
ature to rise and shift slightly toward the inlet of the reactor. This is due to 
the inhibiting action of hydrazine on the rate of ammonia dissociation. Temper- 
ature profiles are shown in Fig. 21 for mass flow rates of 1.5, 3*0, and 6.0 lb/ft~- 
sec. Increasing flow rate causes the peak temperature to rise and shift slightly 
away from the reactor inlet . The effect of changing the catalyst bed configuration 
on the temperature distribution is shown in Fig. 22. Temperature profiles are 
plotted for beds packed with all 25-30 mesh particles, all l/8 in. x 1/8 in . cylin- 
drical pellets, and the standard non-uniform particle size distribution. It is 

apparent that the larger particles slow down the rates of the catalytic decomposition 
reactions. The effects of distributed injectors on temperature profiles are illus- 
trated in Fig. 23. Temperature distributions are plotted for the reference case, 
for the case in which 2/3 of the hydrazine is injected at the inlet and the remaining 
1/3 is injected uniformly over the first l/2 in. of the reactor, and for the case 

in which 1/3 of the hydrazine is injected at the inlet and the remaining 2/3 is 

injected uniformly over the first l/2 in. of the reactor. In these cases the bed 
was taken to be packed with all 25-30 mesh particles . 
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bteady-state temperature and reactant concentration distributions 


a 1 cr 


be obtained using the transient program by running the calculations out to long 
(essentially infinite) times. A comparison of the steady-state axial temperature 
profiles computed from the steady-state and transient programs is shown in Fig. 2k 
for reactor operating conditions of 6= 3-12 lb/ft 2 - sec^ p(feed pressure) = 

V79.5 psia, and T F = 509 deg R and for the standard bed configuration. The rate 
the catalytic decomposition of ammonia, a NHs , was taken as 10^ 
(sec) (see section on Kinetics Information). A comparison of the 
associated axial profiles of mole-fraction of ammonia computed from the steady- 
state and transient models is shown in Fig. 25. 


constant foi 
(lb/it 3) °' t 


Transient Results 


A series of calculations was made using the transient program for the same 
cases cited previously in order to examine the effectiveness of the transient model 
of the reactor system. The results of these calculations are shown in Figs. 2 6 

through 44. In each of these cases the catalyst bed configuration was taken as 
the standard configuration, the feed temperature varied slightly between approx- 
imately 510 and 530 deg R, and the initial bed temperature was taken equal to the 
feed temperature. Temperature distributions are plotted for various times up to 
steady-state in Fig. 2 6 for a case in which the reactor operating conditions were 
taken as G = 3-12 lb/ft 2 - sec and p(feed pressure) = 479-5 psia. The vaporiza- 
tion time noted in Fig. 26 is defined as the time required, under steady-state con- 
ditions, for liquid hydrazine to flow from the reactor inlet to the interface 
between the liquid-vapor and vapor regions. Temperatures are cross-plotted as a 
function of time at three axial positions in Fig. 27* Also shown in Fig. 27 are 
temperature measurements obtained by Rocket Research Corporation (Ref. l6) during 
the course of engine firings under the same operating conditions. The calculated 
mole-fraction profiles corresponding to the temperature profiles shown in Fig. Lx 
are illustrated in Figs. 28, 29, 30 and 31 for hydrazine, ammonia, nitrogen, and 
hydrogen, respectively. The chamber exit pressure for this case was calculated n, 
a function of time and plotted in Fig. 32. It is interesting to note that the cham- 
ber exit pressure rises sharply, initially, as hydrazine and the products of hyd: 
zine decomposition in the upstream end of the reaction chamber flow through the 
reactor, drops as hydrazine decomposes on the catalyst surfaces further downstream 
in the reactor and as the products of hydrazine decomposition diffuse into the 
catalyst particles, and rises again once the concentrations of the products of 
hydrazine decomposition equilibrate between the vapor phase and the catalyst par- 
ticles in the downstream end of the reactor. The results of similar calculations 
made for G = 2.43 lb/ft 2 -sec and P (feed pressure) = 1040 psia are shown, except for 
the concentration profiles, in Figs. 33 through 35 > those made for G= 6.29 
lb/ft 2 -sec and P = 974 psia are shown in Figs. 36 through 38; those made for 
G = 1.52 lb/ft 2 -sec and P = 217.9 psia are shown in Figs. 39 through 4l, and 
those made for G = 1.51 lb/ft 2 -sec and P = 111.4 psia are shown in Figs. 42 
through 44. 
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i t n - xi-c> . _ j.. . r» ■_ _ t_ _ j_ • _ a • j_ • j i . i . — 

xne eiieci/S ox various reauuur operating conuiuions on me variation oi gas 
temperature at the chamber exit with time can be better illustrated by comparing, 
with a reference case, the results of calculations made in turn for various chamber 
pressures, mass flow rates, catalyst bed configurations, axial injection profiles, 
and initial bed temperatures. A reference case was chosen in which the operating 
conditions were taken as G = 3-0 lb/ft 2 - sec, P (feed pressure) = 100 psia, 

Tf r 530 deg R and Tp (initial) = 530 deg R, and the catalyst bed configuration was 
taken as the "standard bed configuration." The bed length was taken as 0.25 ft. 

In Fig. 45 exit gas temperatures are plotted as a function of time for feed pres- 
sures of 100 and 1000 psia with all other conditions taken as those of the reference 
case. Very little effect of pressure on transient response is noted. The marked 
effect of mass flow rate on transient response time is shown in Fig. 46 where exit 
gas temperatures are plotted versus time for mass flow rates of 1.5, 3.0, and 6.0 
lb/ft 2 - sec. The effect of changing the catalyst bed configuration on the temper- 
ature variation with time is illustrated in Fig. 47. Exit gas temperatures are 
plotted as a function of time for beds packed with all 25-30 mesh particles, all 
1/8 in. x 1/8 in. cylindrical pellets, and the standard non-uniform particle size 
distribution. The effects of distributed injectors on temperature profiles is 
shown in Fig. 48. Exit gas temperatures are plotted as functions of time for a 
constant mass flow rate and for the case in which l/3 of the hydrazine is injected 
at the inlet and the remaining 2/3 is injected uniformly over the first l/2 in. 

of the reactor. In these cases the bed was taken to be packed with all 25-30 mesh 
particles and the initial bed temperature was taken as 1200 deg R in order to 
insure vaporization of the cold hydrazine injected downstream of the inlet . 

Finally, the effects of varying the initial bed temperatures are illustrated 
in Fig. 49. Exit gas temperatures are plotted versus time in Fig. 49 for 
initial bed temperatures of 530, 1200, and 1500 deg R. 
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a Radius of spherical particle, ft 

A p Total external surface of catalyst particle per unit volume, of bed, ft" 1 

Cj Reactant concentration in interstitial fluid, Tb/ft^ 

c p Reactant concentration in gas phase within the porous particle, Tb/ft^ 

C p * Equals C p - (c p ) s , lb/ft^ 

Cf Specific heat of fluid in the interstitial phase, Btu/lb - deg R 

Cf Average specific heat of fluid in the interstitial phase, Btu/lb - deg R 

C s Specific heat of catalyst particle, Btu/lb - deg R 

D j Diffusion coefficient of reactant gas in the interstitial fluid^ ft^/sec 

Dp Diffusion coefficient of reactant gas in the porous particle ^ ft^/sec 

F Rate of feed of hydrazine from distributed injectors into the system, 

lb/ft3- sec 

g Dependent variable in Eq. (37) 

G Mass flow rate, lb/ft^- sec 

n Enthalpy, Btu/lb 

A Heat transfer coefficient, Btu/ft^- sec - deg R 

H Heat of reaction (negative for exothermic reaction), Btu/lb 

k c Mass transfer coefficient, ft/sec 

k 0 Reaction rate constant, equals a e Y , sec -1 

Kp Thermal conductivity of the porous catalyst particle, Btu/ft-sec-deg R 

M Molecular weight, lb/lb mole 

M Average molecular weight, lb/lb mole 

n Order of decomposition reaction 

P Chamber pressure, psia 


21 



LIST OF STMBOIS (cont'd) 


A — J- -? • 'D4--.1 /"I V. 

i-iu o_l v ttoxuil ciicx j u o uy a. u iiioxc 

Rate of (heterogeneous) chemical reaction on the catalyst surfaces, 
lb/ft 3- sec 

Rate of (homogeneous) chemical reaction in the interstitial phase, 
lb/ft 3 - sec 

Independent variable in Eq. (37) 

Gas constant, equals 10-73 psia - ft 3 /lb mole - deg R 
Time, sec 

Actual time minus time required, under steady-state conditions, for 
liquid hydrazine to flow from the reactor inlet to the interface between 
the liquid-vapor and vapor regions, sec 

Temperature, deg R 

Mathematical function (defined in Appendix i) 

Mathematical function (defined in Appendix I) 

Weight fraction of reactant in interstitial phase 

Radial distance from the center of the spherical catalyst particle, ft 
Axial distance, ft 

Preexponential factor in rate equation, or 
Constant in Eq. ( 37 ) 

Intraparticle void fraction 

r , ~ ( c p)s H D P 

Kp [ I pj s 

Constant in Eq. (37) 

Equals Q/R(T p ) s 
Interparticle void fraction 

Viscosity of interstitial fluid, lb/ft - sec 
Density of interstitial fluid, lb/ft 3 

Bulk density of catalyst particle, lb/ft 3 
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LIST OF SYMBOLS ( cont * d ) 


«P 


AAo -? /~%r^ -Pi-tv> o*f* t r>n / noH *iri 

A-J4-4, ^XXV^lliCAUJ.'^Ci-L. XUiiUVXVIX ^ LLw -L -LHC wi. JLI1 


A r\r\Qn rl -i -v 

■» X J.V-4— 


T 

■*■ J 


Subscripts 


F Refers to feed 

i Refers to interstitial phase } or 

Time index in transient model 
j Axial distance index in transient model 

k Index in Eqs . ( 38 ) and ( 39 ) 

P Refers to gas within the porous catalyst particle 

s Refers to surface of catalyst particle 


Superscripts 

J Refers to chemical species 

L Refers to liquid at vaporization temperature 

V Refers to vapor at vaporization temperature 
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APPENDIX 1 


DEP.r/ATION OF INTEGRAL EQUATIONS REPRESENTING THE CONCENTRATION 
PROFILES OF REACTANTS WITHIN THE CATALYST PARTICLES 


In this section equations are developed to describe the steady-state concen- 
tration profiles of hydrazine vapor and of ammonia within the catalyst particles. 
The reactant concentration profiles in the porous particles at any axial location 
can be found as solutions to: 


D P V C P ' r het IC P ) = 0 


( 1 - 1 ) 


If the catalyst particles are taken to be "equivalent" spheres of radii a, and if 
concentration Cp* is defined such that = C p - (Cp) $ , Eq. (l-l) can be written as 


j d_ 

x 2 dx 


d c 

_ r 


* \ “1 


tel ' 0 


( 1 - 2 ) 


where x is the radial distance from the center of a sphere. The boundary conditions 
associated with Eq. (1-2) are 


= 0 AT X = 0 , 


d c p 
dx 


= 0 AT x = 0 


Equation (l-2) can be rearranged to get 


dx 


d c p 
dx 


C ,x £ 




Dp 


x, c p 


(1-3) 


(1-4) 


The solution to Eq. (1-4) is most easily obtained by converting it into a Fredholm 
integral equation (see Ref. 17 ) of the form 


C p *(x) - 


x 2 [u(x) v\x) - t/uivixjj o 


/ a G(x,C) ^>(C,Cp *)di ( I -5) 
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• where u (x) is a solution of 


dx \ d* / 

subject to the condition that 


= 0 

x = o 


dc£ 

dx 


du * 
dx Cp 


and v (x ) is a solution of 




subject to the condition that 


d C 


dx 


d\j * 

~d7 C P 


= o 


u x= a 


The Green's function, G (x,£) is given by 


f u(£) v (x) FOR 0 < ( < x 
\ulx) v(£) FOR x < £ < a 


( 1 - 6 ) 


(1-7) 


(1-8) 


(1-9) 


( 1 - 10 ) 


The function u (x) can be determined by first integrating Eq. (1-6) to get 


U 


+ B, 


( 1 - 11 ) 


Applying Eq. (1-7) together with the first of boundary conditions (1-3) to 
Eq. (1-11), it is found that A - 0 and 


u = B, 


( 1 - 12 ) 
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The function v(x) can be determined in a similar manner by first integrating 
Eq. (1-8) to get 


v = 



(1-13) 


and then applying Eq. (1-9) and the second of boundary conditions (1-3) to Eq. 
(1-13) to get 



(1-14) 


Equations (1-10), (l-12), and (l-l4) can now be combined to get 


GU,£ ) 


f A 0 B 

■J_ _ i ' 

) 2 1 
\ 

.a ~ x . 

A„B, 

'j i_" 

^ 2 1 

a i- 


FOR 

FOR 


0 < £ < x 
x < C < a 


In addition, 


u 


(x) V (x) - 


u'(x) 


v(x) 




A 2 B , 


(1-15) 


(1-16) 


Equations (l-15) and (l-l6) can now be substituted into Eq. (l~5) to get 


C 











(1-17) 


or 


Cp (x) 





r he» (C P 


dt 




r het 


di 


(1-18) 


Equation (l-l8) is an implicit integral equation which can be solved numerically 
to determine the concentration at any point in a porous particle in terms of (Cp) s , 
the concentration at the surface of the particle. The gradients of Cp, evaluated 
at the particle surfaces at given axial locations, can then be obtained from the 
slopes of curves of Cp versus x at X = Q. 
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NORMALIZED RADIAL DISTANCE WITHIN PARTICLE, x/Q 



EFFECT OF REACTION RATE CONSTANT 
ON TEMPERATURE PROFILE WITHIN CATALYST PARTICLE 



NORMALIZED RADIAL DISTANCE WITHIN PARTICLE, x/Q 




VARIATION OF CONCENTRATION GRADIENT AT PARTICLE SURFACE 

WITH REACTION RATE CONSTANT 



REACTION RATE CONSTANT 
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STEADY-STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS RATE CONSTANTS FOR THE CATALYTIC 
DECOMPOSITION OF AMMONIA 
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MOLE- FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 6 


STEADY-STATE AXIAL PROFILES 
OF MOLE-FRACTIONS OF REACTANTS 


P = 479.5 PSIA 
G = 3.12 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




FRACTIONAL AMMONIA DISSOCIATION 


FIG. 7 




TEMPERATURE IN INTERSTITIAL PHASE, T; - DEG 
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FIG. 8 


STEADY-STATE AXIAL TEMPERATURE PROFILE 


P = 1042 PSIA 
G = 2.43 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




MOLE - FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 9 


STEADY -STATE AXIAL 
OF MOLE -FRACTIONS OF 


PROFILES 

REACTANTS 


P = 1042 PSIA 
6 = 2.43 LB/ FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




FRACTIONAL AMMONIA DISSOCIATION 
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FIG. 10 


STEADY-STATE AXIAL PROFILE 
OF FRACTIONAL AMMONIA DISSOCIATION 


P = 1042 PSIA 
G = 2.43 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, T, - DEG R 
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FIG. II 
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AXIAL TEMPERATURE PROFILE 


P = 974 PS! A 
G = 6.29 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION {SEE TEXT) 




MOLE-FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 12 


STEADY-STATE AXIAL PROFILES 
OF MOLE- FRACTIONS OF REACTANTS 


P = 974 PSIA 
G = 6.29 LB/FT 2 - SEC 

STANOARO BED CONFIGURATION (SEE TEXT) 
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FRACTIONAL AMMONIA DISSOCIATION 
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FIG. 13 


STEADY-STATE AXIAL PROFILE 
OF FRACTIONAL AMMONIA DISSOCIATION 


P = 974 PS I A 
G = 6.29 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, T; - DEG 
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FIG. 14 


STEADY -STATE AXIAL TEMPERATURE PROFILE 


P = 217.9 PSIA 
G = 1.52 LB/ FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




MOLE - FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 15 


STEADY-STATE AXIAL PROFILES 
OF MOLE-FRACTIONS OF REACTANTS 


P = 217.9 PSIA 
G = 1.52 LB /FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 




FRACTIONAL AMMONIA DISSOCIATION 
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FIG. 16 


STEADY -STATE AXIAL PROFILE 
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P = 2179 PSIA 
G = I .52 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 








MOLE-FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 18 


STEADY-STATE AXIAL PROFILES 
OF MOLE- FRACTIONS OF REACTANTS 


P = III. 4 PSIA 
G = 1.51 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




FRACTIONAL AMMONIA DISSOCIATION 
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FIG. 19 


STEADY-STATE AXiAL PROFILE 
OF FRACTIONAL AMMONIA DISSOCIATION 


P = 1 1 1.4 PSIA 
G = 1.51 LB /FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 


F9I046I-I2 F| G. 20 

STEADY-STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS CHAMBER PRESSURES 

G = 3.0 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 





TEMPERATURE IN INTERST ITI AL PHASE, T j - DEG 
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FIG. 22 


STEADY-STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS CATALYTIC BED CONFIGURATIONS 

P = IOO PSIA 
G = 3.0 LB/FT 2 - SEC 




TEMPERATURE IN INTERSTITIAL PHASE, T: - DEG 
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F9 1 046 1--12 FIG. 24 

COMPARISON OF STEADY-STATE AXIAL TEMPERATURE PROFILES 
COMPUTED FROM STEADY -STATE AND TRANSIENT MODELS 


P = 479.5 PSIA 
G = 3.12 LB /FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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FIG. 25 


COMPARISON OF STEADY-STATE AXIAL PROFILES 
OF MOLE - FRACTION OF AMMONIA 
COMPUTED FROM STEADY- STATE AND TRANSIENT MODELS 

P = 479.5 PSIA 
G =3.12 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 
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TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 
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MG. 27 


VARIATION OF TEMPERATURE WITH TIME 
AT SEVERAL AXIAL POSITIONS 


P (FEED PRESSURE) = 479.5 PS I A G = 3.12 LB/FT*-SEC STANDARD BED CONFIGURATION (SEE TEXT) 
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MOLE- FRACTION OF HYDRAZINE IN INTERSTITIAL PHASE 
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TRANSIENT AXIAL PROFILES 
OF MOLE -FRACTION OF HYDRAZINE 


FIG. 28 


P ( FEED PRESSURE) = 479.5 PSIA 
G = 3.12 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 





MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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TRANSIENT AXIAL PROFILES 
OF MOLE-FRACTION OF AMMONIA 


FIG. 29 


P(FttL) pRESSUnE) = 473.5 PS I A 
G = 3.12 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




MOLE - FRACTION OF NITROGEN IN INTERSTITIAL PHASE 
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FIG. 30 


TRANSIENT axial profiles 
OF MOLE -FRACTION OF NITROGEN 


P ( FEED PRESSURE) - 479.5 PSIA 
G = 3.12 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 





MOLE- FRACTION OF HYDROGEN IN INTERSTITIAL PHASE 
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TRANSIENT AXIAL PROFILES 
OF MOLE-FRACTION OF HYDROGEN 


MG. 31 


P (httL) PRESSURE) = 479.5 PSIA 
G = 3.12 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 
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FIG. 32 


VARIATION OF CHAMBER EXIT PRESSURE WITH TIME 


P ( FEED PRESSURE) = 479.5 PSIA 
G =3.12 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 



TEMPERATURE IN INTERSTITIAL PHASE, T, - DEG R 
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FIG. 33 


TRANSIENT AXIAL TEMPERATURE PROFILES 

P (FEED PRESSURE) : !042 PSIA 
G = 2.43 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 





TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG P 
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FIG. 34 


VARIATION OF TEMPERATURE WITH TIME 
AT SEVERAL AXIAL POSITIONS 


P (FEED PRESSURE) = 1042 PSIA G = 2.43 LB/FT 2 - SEC STANDARD BED CONFIGURATION (SEE TEXT) 
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CHAMBER EXIT PRESSURE - PSIA 
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FIG. 35 


VARIATION OF CHAMBER EXIT PRESSURE WITH TIME 


►'rctaaursc ; - ium^ rain 
G = 2.43 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 



TIME, t - SEC 




TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 
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FIG. 36 


TRANSIENT AXIAL TEMPERATURE PROFILES 

P ( FEED PRESSURE) = 974 PSIA 
G = 6,29 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 
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FIG. 


VARIATION OF TEMPERATURE WITH TIME 
AT SEVERAL AXIAL POSITIONS 


P (FEED PRESSURE) = 974 PSIA G = 6.29 LB/FT 2 - SEC STANDARD BED CONFIGURATION (SEE TEXT) 
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CHAMBER EXIT PRESSURE - PSIA 
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MG . 38 


VARIATION OF CHAMBER EXIT PRESSURE WITH TIME 


P ( FEED PRESSURE) = 974 PSIA 
G = 6.29 LB/FT 2 -SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, T; - DEG 
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MG . 39 


I KANSiENT 


AA I ML 
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P ( FEED PRESSURE) - 217.9 PSIA 
G = 1.52 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 
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VARIATION OF TEMPERATURE WITH TIME 
AT SEVERAL AXIAL POSITIONS 


FIG. 40 


P i TEED PRESSURE) = 217.9 PSIA G = 1.52 LB/FT 2 - SEC STANDARD BED CONFIGURATION (SEE TEXT) 
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PRESSURE WITH TIME 


• «/ < 


G = 1.52 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 
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FIG. 42 
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TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 
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FIG. 4 3 


VARIATION OF TEMPERATURE WITH TIME 
AT SEVERAL AXIAL POSITIONS 


P (FEED PRESSURE) = III. 4 PSIA G = 1.51 LB/FT 2 - SEC STANDARD BED CONFIGURATION (SEE TEXT) 
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CHAMBER EXIT PRESSURE - PSIA 
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MG. 44 


VARIATION OF CHAMBER EXIT PRESSURE WITH TIME 


P ( FEED PRESSURE) = 111.4 PSIA 
G = 1.51 LB/FT 2 - SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 




EXIT GAS TEMPERATURE - DEG R 




EXIT GAS TEMPERATURE - DEG 
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FIG. 46 


VARIATION OF EXIT GAS TEMPERATURE WITH TIME 
FOR VARIOUS MASS FLOW RATES 


P ( FEED PRESSURE) = IOO PSIA 
STANDARD BED CONFIGURATION (SEE TEXT) 




EXIT GAS TEMPERATURE - DEG R 
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VARIATION OF EXIT GAS TEMPERATURE WITH TIME 
FOR VARIOUS CATALYTIC BED CONFIGURATIONS 

P ( FEED PRESSURE) = IOO PSIA 
G = 3.0 LB/FT 2 - SEC 





EXIT GAS TEMPERATURE - DEG R 
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VARIATION OF EXIT GAS TEMPERATURE WITH TIME 
FOR VARIOUS HYDRAZINE AXIAL INJECTION PROFILES 


P ( FEED PRESSURE) = 100 PSIA 
BED CONFIGURATION: ALL 25 - 30 MESH GRANULES 
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VARIATION OF EXIT GAS TEMPERATURE WITH TIME 
FOR VARIOUS INITIAL BED TEMPERATURES 


P(FEED PRESSURE) = IOO PSlA 
G * 3.0 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 
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